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Dissipative quantum systems are sometimes phenomenologically described in terms of a non- 
hcrmitian hamiltonian H , with different left and right eigenvectors forming a bi-orthogonal basis. 
It is shown that the dynamics of waves in open systems can be cast exactly into this form, thus 
providing a well-founded realization of the phenomenological description and at the same time 
placing these open systems into a well-known framework. The formalism leads to a generalization 
of norms and inner products for open systems, which in contrast to earlier works is finite without 
the need for regularization. The inner product allows transcription of much of the formalism for 
conservative systems, including perturbation theory and second-quantization. 
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Introduction 

Dissipative systems can be discussed in many ways. 
The fundamental approach recognizes that energy flows 
from the system S to a bath B, whose degrees of freedom 
are then eliminated from the path integral or equations 
of motion jjj . While rigorous, this approach is inevitably 
complicated, and often leads to integro-differential equa- 
tions for time evolution. An alternate phenomenological 
approach postulates a non-hermitian hamiltonian (NHH) 
H, whose left and right eigenvectors form a bi-orthogonal 
basis (BB) §-§]. These NHHs with discrete BBs can 
sometimes be obtained from a full quantum theory, but 
usually under some approximations ■ 

This Letter discusses a class of models of waves in open 
systems. These are scalar fields (j>{x, t) in 1 d, described 
by the wave equation. Outgoing wave boundary condi- 
tions cause the system to be dissipative. We show that 
these open systems are exactly described by an NHH with 
a BB formed by the resonances or quasinormal modes 
(QNMs). This connection on the one hand provides the 
phenomenological approach with a realization which has 
an impeccable pedigree rigorously traceable to the fun- 
damental approach, and on the other hand places earlier 
work on such open systems into a familiar framework. 
A generalized inner product emerges; in contrast to pre- 
vious works, it is finite and requires no regularization. 
Under the generalized inner product, the hamiltonian H 
is symmetric, which opens the way to a clean formulation 
of perturbation theory and second-quantization in terms 
of the QNMs of the system. 
Waves in Open Systems 

We consider waves in I d described by 
[p(x)df — <9 2 ] 4>{x,t) = on the half line [0, oo), with 
4>(x = 0, t) = and 4>(x,t) approaching zero rapidly as 
x — > oo P). Let the system S be the "cavity" / = [0,a], 
and the bath B be (a, oo), where p{x) = I. Energy is 



exchanged between S and B only through the bound- 
ary x — a. We impose the outgoing wave condition 
dt4>{x,t) — ~d x (j){x,t) for x > a. 

This mathematical model is relevant for many phys- 
ical systems: the vibrations of a string with mass den- 
sity p the scalar model of EM in an optical cavity 
(the node at x = is a totally reflecting mirror, and 
a partially transmitting mirror at x = a can be mod- 
eled by p(x) = MS(x — a)) |pTj| ; or gravitational radi- 
ation from a star with radius a fl^| . The wave equa- 
tion can be mapped to the Klein-Gordon equation with 
a potential V(x) [jl3), which is relevant for gravitational 
waves Jjl|] ; here <fi is the perturbation about the spherical 
background metric of a star, a; is a radial coordinate re- 
lated to the circumferential radius r, and V describes the 
wave scattering by the background metric. Gravitational 
waves carrying the signature of the QNMs of black holes 
may soon be observed by new detectors such as LIGO 
and VIRGO @. 

For the "cavity" / = [0, a], the outgoing condition 
is imposed at x = a + only. The QNMs are factor- 
ized solutions on /: 4>(x,t) — f n (x)e~ ZUJnt , with [<9 2 + 
p(x)uj^ l ] f n {x) = 0. These are observed in the frequency 
domain as resonances of finite width (e.g., the EM spec- 
trum seen outside an optical cavity) or in the time do- 
main as damped oscillations (e.g., the numerically sim- 
ulated gravitational wave signal from the vicinity of a 
black hole). It would obviously be interesting to be able 
to describe these QNMs in a manner parallel to the nor- 
mal modes (NM) of a conservative system. 

These QNMs form a complete set on I if (a) p(x) has 
a discontinuity at x — a to provide a natural demar- 
cation of the "cavity", and (b) p(x) = 1 for x > a, 
so that outgoing waves are not scattered back into the 
system |16| . Under these conditions, one can expand 
4>{x,t) — J2n a nfn{x)e~ luJnt for x £ I and t > 0, thus 
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allowing an exact description of the system in terms of 
discrete variables (modes spaced by Alo ~ rather 
than a continuum. Nevertheless, the analogy with con- 
servative systems is still not apparent: Is there a natural 
inner product (with which to do projections and thus to 
prove the uniqueness of expansions)? Is there a norm to 
scale wavefunctions (noting that f n diverges at spatial in- 
finity)? Can perturbation theory be formulated (noting 
that the usual proofs require an inner product to define 
orthogonality)? Can the theory be second-quantized? 
This Letter shows that all these questions have natural 
answers in the language of a BB. 

Phenomenological non-Hermitian Hamiltonians 
and Bi-orthogonal Bases 

Though not rigorously founded upon a genuine quan- 
tum theory, NHHs with BBs are nevertheless well de- 
veloped as a postulatory system ^,|). Consider a 
space W on which is defined a non-hcrmitian opera- 
tor H and a conjugate linear duality transformation 
D: D (a|$)+ /3|*)) = a* D |$) + (3* D |#), such that 
DH = H^D fjj. The BB consists of the two set of 
eigenvectors \F n ) G W and \G n ) = D \F n ) G W = D(W) 
satisfying H\F n ) = w n \F n ), W \G n ) = u* \G n ), where 
the two eigenvalues are related by duality. By project- 
ing the eigenvalue equations on (G n \ and \F n ), it follows 
easily that (G n \F m ) = 0, for m ^ n. 

It is usually assumed that these eigenstates are com- 
plete, so that any vector can be expanded as |$) = 
Y,n a n \ F n), wi th a n = (G n \$) / {G n \F n ) , leading imme- 
diately to the resolution of the identity and of the time- 
evolution operator 
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which in principle solves all the dynamics |TJ] . 
Bi-orthogonal Basis for the Wave Equation 

BBs are widely used in many disciplines, for example in 
the theory of wavelets [|l9| and to describe excited molec- 
ular systems [p|j20|. The left and right eigenvectors of 
the Maxwell operator are typically used to represent the 
Green's function for EM fields in open cavities [^l|-p3|, 
or to evaluate Fox-Li states p4|. Here we seek a paral- 
lel with quantum mechanics, similar to earlier works for 
generalized oscillators (25) and the classical wave equa- 
tion (without dissipation due to leakage) . The prob- 
lem at hand, where there is dissipation due to outgoing 
waves, was formulated in this manner recently p7[ , and 
is briefly sketched below, especially as it relates to the 
BB. 

It is natural to introduce the conjugate momentum cf> — 
p{x)dt4>, and the two-component vector |$) = (</>, ^) T . 
In terms of this, the dynamics can be cast into the 
Schrodinger equation with the NHH 
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The identification <f> = pdt4> follows from the evolution 
equation |28|. 

The natural definition of an inner product between 
|*) = (i/>,-0) T and |$) = (0, cf>) T on [0, oo) is 



(*|4> 



dx 
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However, on account of the assumed asymptotic behav- 
ior, the integral is convergent. 

For outgoing waves, we consider only the space U of 
such vectors |<E>) defined on [0, oo) which satisfy the out- 
going condition <f> — ~(j)' for x > a. The bath variables 
are eliminated simply but exactly by projecting to the 
space W of vectors |$) defined on / and which satisfy 
4> = —ft at x = a + . The QNMs are right-eigenvectors 
of F: \F n ) = (/„,/„) T = (fn,-iu n pf n ) T . The duality 
transformation is Z3(0i,02) T = — £($j> ( / > i) T - 

For open systems, a crucial concept is the inner prod- 
uct between one vector and the dual of another, to which 
we give a compact notation: 



ip4>) dx 
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which is linear in both vectors, and cross-multiplies the 
two components, properties to be emphasized below. 
This bilinear map plays the role of the inner product 
for conservative systems. 

Our notation does not distinguish between functions 
(say |$)) defined on [0, oo) and their restrictions to /; 
the former are in U and the latter are in W, with the as- 
sociation between them being many-to-one. As written 
in (||), the inner product involves the wavefunctions out- 
side /, i.e., it appears to be defined on U rather than W. 
However, one can completely eliminate the bath degrees 
of freedom: because of the outgoing conditions, the inte- 
grand on (a, oo) reduces to a total derivative, and (ph can 
be written purely in terms of the inside variables p7| : 




V>0 dx + i)(a + )(f>{a + ) ) (6) 



The surface term is the only remnant of the outside. 
Thus, (|^) can be regarded as a bilinear map (or loosely 
an inner product) defined on W p7| . The somewhat 
peculiar structure (e.g., the cross- multiplication between 
the two components and the appearance of the surface 
term) is now seen to arise naturally from (^) upon the 
introduction of the duality transformation. In the limit 
where the escape of the waves is small, the generalized 
norm of an eigenvector [F n , F n ) reduces to 2w n times the 
conventional norm; this is the reason for choosing the 
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phase convention for D. The ability to normalize QNM 
wavefunctions is nontrivial, since /„ diverges at spatial 
infinity, and a naive expression such as / °° \f n \ 2 dx would 
not be appropriate. 

The diagonal version (<E>, <&) for the special case of 
QNMs was first introduced by Zeldovich |29| in a form 
that involved (a) <fi outside / (so that it is defined on U 
rather than W) and (b) regularization of the divergent in- 
tegral rather than a surface term; it was later re-cast into 
the form and generalized to 3 d and EM fields [ 30 1 . 
The off-diagonal form ('I',^) was later introduced [27|. 
Here, by relating the discussion to bi-orthogonal states 
and the duality transformation, it is seen that these con- 
cepts emerge naturally, including the specific form of (j^) . 

An inner product equivalent to @) has also been dis- 
cussed extensively from other perspectives 31 3^]. In 
these works, the inner product is defined on [0, oo) rather 
than a finite interval, with the consequent divergence 
(e.g., for the inner product between two QNMs each 
growing exponentially at infinity) handled either (a) by 
a regulating factor exp(— ex 2 ), e — > + , (b) analytic con- 
tinuation in the wavenumber k, or (c) complex rotation 
in the coordinate x. Each of these procedures has its 
limitations; in contrast, (JsJ) makes no reference to the 
outside or bath, and is computationally convenient and 
manifestly finite. 

Under this bilinear map, H is symmetric: (*f?,H$>) = 
(^ji?^), which follows very simply from DH — H^D. 
This key property is analogous to the hermiticity of H 
for conservative systems. It is nontrivial, in that surface 
terms that arise in the integration by parts are exactly 
compensated by the surface terms in (^). This symmetry 
property leads, in the usual way, to the orthogonality of 
non-degenerate eigenfunctions. 

The completeness relation (0) is a dyadic equation. Its 
(1, 2) and (1, 1) components lead to the sum rules |34] 
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E \fn{x)fn{y)p{x) = iS(x - y) 



(7) 



for x,y £ I, which have been derived and discussed ex- 
tensively J27|. 

The completeness and orthogonality relationships es- 
tablish the QNMs as a BB, and moreover allow the time 
evolution to be solved as \&(x,t)) = J2 n °n e ~' w °' l-^n), 
where a n — (G n \^>(x 1 0))/(2Lu n ). This is a discrete and 
exact representation of the dynamics, even though / is 
open to an infinite universe with a continuum of states. 
Completeness is not proved in most other applications of 
NHHs to physical systems. 
Perturbation theory 

These notions allow much of the standard formalism 
in quantum mechanics to be carried over. As one exam- 
ple consider time-independent perturbation theory. Let 



po(x)^ 1 be changed to p{x)^ 1 — po(x)^ 1 [1 + pV(x)], 
where \fi\ <C 1 V(x) has support in /. Then the per- 
turbation to the eigenvalues and eigenfunctions can be 
written in the standard Rayleigh-Schrodinger form, in 
terms of a discrete series p7| . These formulas, though 
superficially identical with textbook formulas for conser- 
vative systems, are nontrivial in two ways. First, the 
perturbative formulas apply to complex eigenvalues. Sec- 
ond, the use of resonances implies that there is no "back- 
ground" , and expressing the corrections in terms of dis- 
crete modes also means that the small parameter of ex- 
pansion is /x/|Aw| ~ fia/n, which would not have been 
apparent in terms of the states of the continuum. 

The derivation of these results simply follows the con- 
servative case (everywhere replacing inner products by 
the bilinear map (\&, $)), and need not be repeated. 
Discussion 

We have established an exact correspondence between 
phenomenological NHHs and waves in a class of open 
systems. This relationship provides a well-founded real- 
ization of NHHs. Because we start with a hamiltonian 
system and remove the bath degrees of freedom with- 
out approximations, these open systems can be second- 
quantized |3jJ. In other words, one can discuss photons 
in open cavities using BBs, which makes this class of 
examples unique and interesting. The relationship also 
places these open systems into a well-known and conve- 
nient framework. Thus, the linear space structure, or- 
thogonality and completeness can all be derived natu- 
rally, by transcribing usual derivations for conservative 
systems and everywhere replacing the inner product by 

The formalism discussed here also applies to the Klein- 
Gordon equation with a potential V(x) which ap- 
plies, among other things, to linearized gravitational 
waves propagating away from a black hole. The first- 
order perturbation result for the QNM frequencies has 
been used to understand the shifts in the gravitational 
wave frequencies when a black hole is surrounded by an 
accretion shell |33|] . 

The wave equation discussed here may be regarded as 
a physical realization of BBs for open systems. Many 
other inequivalent realizations arise when one considers 
outgoing waves in a spherically symmetric 3-d system; 
each angular momentum I leads to realizations in which 
the surface terms in the inner product involves I radial 
derivatives |36| ]. 

However, the entire formalism refers to systems de- 
scribed by second-order differential equations, so that 
two sets of initial data, namely (f> and (f>, are required, 
and the outgoing condition is expressed as a constraint 
between them. The formalism does not apply in its 
entirety to systems described by first-order differential 
equations, e.g., a-decays described by the Schrodinger 
equation with Gamow boundary condition. In any event, 
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the Schrodinger equation formally gives unbounded sig- 
nal speeds and does not possess outgoing and incom- 
ing sectors related by time reversal; thus the concept of 
outgoing waves is actually quite different. Nevertheless, 
if one is interested only in frequency domain problems, 
e.g., eigenvalue problems and time-independent pertur- 
bation theory, then the formalism survives even in this 
case. This is most easily appreciated by starting with the 
Klein-Gordon equation and simply relabelling lo 2 h- > uj. 

Using ( V E',$) rather than the equivalent form (D^l^) 
allows all reference to D to be avoided. However, ("J, $) 
is a bilinear map (rather than being linear in the ket 
and conjugate linear in the bra). This property is quite 
general, since D is conjugate linear. But in most applica- 
tions of the inner product (e.g., for projections), it does 
not matter whether the map is linear or conjugate linear 
in the bra; this is why results from conservative systems 
can be carried over. The only property that is lost is the 
positivity of (<&,$), which is unsurprising for a dissipa- 
tive system. Thus it is useful to think of the states of 
quantum dissipative systems as vectors in a linear space 
W endowed with such a bilinear map, which is the gener- 
alization of the notion of a Hubert space. Time-evolution 
is then generated by an operator H which is symmetric. 

The open systems described here are genuinely dis- 
sipative, with Im io n < 0. This contrasts with some 
models with NHHs which are nevertheless conservative 
P^p6| . For infinite-dimensional NHH models, complete- 
ness of the BB is usually assumed, but difficult to prove. 
Through these wave systems, we have provided explicit 
examples where completeness can be proved (if the dis- 
continuity and "no tail" conditions are met), as well as 
examples where the basis is not complete (if these con- 
ditions are not met). These should also be useful in fur- 
thering understanding of NHH models. 
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